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The ring configurations for classical two-dimensional (2D) atoms are calculated within the Thom- 
son model and compared with the results from 'exact' numerical simulations. The influence of the 
functional form of the confinement potential and the repulsive interaction potential between the 
particles on the configurations is investigated. We also give exact results on those eigenmodes of 
the system whose frequency does not depend on the number of particles in the system. 
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I. INTRODUCTION 



Recently, a detailed investigation of the ground state configurations of a finite two-dimensional (2D) classical system 
of charged particles, confined by a parabolic external potential, was made in Ref. [l]. It was found that the particles 
arranged themselves on rings. A study of the spectral properties of this classical system such as the energy spectrum, 
the eigenmodes, and the density of states was made in Ref. |^. In both papers Monte Carlo simulations were used to 
study these classical atoms. 

In the present paper we aim to obtain analytic results for these classical 2D atoms by using a model system- the 
Thomson model. Thomson proposed this classical model in 1904 in order to calculate the structure of the atom.B He 
was unable to obtain analytic results for the case of a real 3D atom and therefore constrained the particles to move 
in a plane. Together with the parabolic confinement potential this is now a model for the classical 2D atoms which 
were studied in Refs. [I] and ||. A crucial ansatz in the Thomson model is that the particles arrange themselves in a 
ring structure. Next, Thomson puts as many particles as possible on a single ring, as allowed by stability arguments, 
and the rest of them are placed in the centre. By successive applying this approach he was able to construct the 
approximate atomic ring configurations. 

Here we apply the Thomson model and compare the obtained configurations with those from the 'exact' numerical 
simulations of Ref. |l| and determine when the Thomson model breaks down. In a second step we generalize the 
Thomson model to: 1) general r n confinement potentials, and 2) different functional forms (1/r™ and logarithmic) 
for the repulsive inter-particle potential. 

This paper is organized as follows. The calculation to obtain the structure and the eigenfrequencies in the Thomson 
model is outlined in Section We compare the results for parabolic confinement and Coulomb repulsion with the 



previous results from 'exact' numerical simulations. Section [II is devoted to the extensions of the Thomson model. 



In Section IV we give some analytical results for the eigenfrequencies of artificial atoms. Our conclusions and our 
results are summarized in Section [v|. 



II. THE THOMSON MODEL 



A. General system 

We study a system of a finite number, TV, of charged particles interacting through a repulsive inter-particle potential 
and moving in two dimensions (2D). A confinement potential keeps the system together. We focus our attention on 
systems described by the Hamiltonian 

i— 1 j>i ' J ' 

where m is the mass of the particle, wq the radial confinement frequency, e the particle charge, e the dielectric constant 
of the medium the particles are moving in, and fl = (xi, yi) the position of the particle with = \fi\. For convenience, 
we will refer to our charged particles as electrons, keeping in mind that they can also be ions with charge e and mass 
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to. We can write the Hamiltonian in a dimensionless form if we express the coordinates, energy, and force in the 
following units 

r = (e 2 /e) 1 - / (™+™') a -V(«+n') ; ( 2 .2a) 

E a = ( e 2 / e ) n /( n + n ') a n ' /(n+n '\ (2.2b) 

Fq = ( e 2/ £ )(»-l)/(™+™') a (n' + l)/(n+n') j p.2c) 

where a = ?muj 2 /\ n . All the results will be given in reduced form, i.e., in dimensionless units. In such units the 
Hamiltonian becomes 

^ = E'-r + E F J - F - (2.3) 

i=l 3 >i 1 * jl 

Thus, the energy is only a function of the number of particles N, the power of the confinement potential n and the 
power of the interaction potential n' . For n — 2 and n' = 1 this system reduces to the one studied in Refs. [I] and ||. 
The Humerical values for the parameters ro and Em for some typical experimental systems like (electrons in quantum 
dotsja electron bubbles on a liquid Helium surface,B and ions trapped in Penning and Paul trapscl were given in Ref.| 
for the case of parabolic confinement and Coulomb repulsion. 

In the Thomson model one obtains the groundstate configurations of this system by making the following assump- 
tion: the N electrons are arranged at equal angular intervals around the circumference of a circle of radius a. Then 
one investigates the stability of this configuration. With other words, small displacements of the electrons out of their 
equilibrium have to remain small in time. Solving the equations of motion in polar coordinates to first order in these 
displacements, we obtain an equation which determines the eigenfrequencies to of the system 

((" + "0 2 »,+i"2+2 S ^ + L k -L + {n- 2)fi 2 - w 2 ) (N - N k - u; 2 ) = (M k - 2flu;) 2 , (2.4) 

where f2 is a constant angular velocity with which the whole system rotates around its centre, I(n') — n', and 



N-l 



S 



N 



I(n') v-^ 1 cos(2kjTT /N) , , , » , , . 

L k = , \ ' V )J ' ' (n'sin 2 (jn/N) + l) , (2.5b 

( 2a )«'+2 ^ S m n +2 (jir/N) K U ' ' h V ' 

H n ') v^ 1 cos(2fc77rA/V) . , . . . . „ . 

N k = - \ ' > \ J ' — '- ((n! + 1 cot 2 (77r/7V) + l) , (2.5c 

(2a)"'+ 2 ^ S in™ (jn/N) U ' W ' ' ' 

M k =n' \ ' . > sin— rrr-^ — —■ (2.5d) 

2a ™'+ 2 ^ N sin n +1 (j7r/iV) V ; 



V ^ , (2.5a) 



where k is an integer between and N — 1. The derivation of Eq. (2.4) is straightforward and proceeds along the 
lines given by ThomsonB for the case of n — 2 and n' = 1. How many frequencies yields this equation? From the 



onset we notice that if we replace k by N — k in Eq. (2.4), the values of lo differ only in sign, and thus results in the 



same frequencies. Consequently all the values of u) can be obtained by taking only k = 0, 1, ... , ( N — l)/2 if N is 



odd, or k = 0, 1, . . . , N/ 2 if N is even. Thus, if N is odd there are (N + l)/2 equations of the type (|2.4| ). For k = 



M k=0 = 0, and Eq. (2.4) reduces to a quadratic equation, which implies that the number of roots of these (N + l)/2 
equations equals 4 x (N + l)/2 — 2 = 2N, i.e. the number of degrees of freedom of the system, as should be. For N 
even, there are N/2 + 1 equations. But M k — for k = and k = N/2, and thus two of these reduce to quadratics. 
Consequently, the number of roots is 4 x (N/2 + 1) — 4 = 2N. Thus in both cases the number of roots is equal to 
2N, the number of deg rees of freedom of the electrons in the plane of motion and therefore all eigenfrequencies are 
obtained from Eq. (fO|). 
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B. Parabolic confinement and Coulomb repulsion 



Le t's c onsider the case of parabolic confinement (n — 2) and Coulomb repulsion between the electrons (n! = 1). 



Eq. (2.4) for the eigenfrequencies becomes 

(J^Sn + L k -L - uA (N - N k - u 2 ) = (M k - 2Quj) 2 . (2.6) 



If we take A = 1 in Eq. (2T) the frequency is expressed in the unit u>' — ljo/^/2. 

Figure [l] shows the eigenfrequencies squared for N — 1,2, 3, 4, 5 and 6 for zero angular velocity, i.e. fl = 0. Notice 
that the lowest non zero eigenfrequency decreases and in fact for N — 6 we find that ui 2 < indicating that the single 
ring configuration is no longer stable. However, it is pos sible to stabilize the system by placing electrons in the centre. 



If we put p electrons in the centre of the ring, Eq. (2.4) is modified into 

f, „ Hn') (n' + 2)I(n')p s , ,\ , 9 , 

((" + " V + i n- + 2 g " + L * ~ L * + J + 2 + ( n ~ 2 )" 2 - - 2 j (N -N k - c 2 ) 

= (M k - 2Quj) 2 . (2.7) 

From this equation we obtain the minimum value for p which is needed in order to make a ring of N out electrons 
stable. For parabolic confinement and Coulomb repulsion between the particles we find the condition 

p > f(k, N out ) ^ a l(Lo~L k )~^ + £ _ML_, (2.8) 

which must be satisfied for every k. We introduce f(N out ) = maxfc g [ .Ar) and p equals the integer which is just larger 
than f(N out ). Figure || shows the function f(N out ). For p larger than one, the p inner electrons in principle cannot 
be in the same point, i.e. in the centre. They will repel each other until they balance the confinement potential. The 
Thomson model assumes, as is approximately the case, that the p electrons around the centre exert the same force 
as resulting from p c har ges placed at the centre. As an example we consider the case of N ou t = 12. For an outer ring 



of 12 electrons, Eq. (2.8) requires that 7 electrons are inside (this result can be read off from Fig. |2|). But 7 electrons 
cannot form a single ring, but will arrange themselves as a ring of 6 with one at the centre. Thus the system of 19 
electrons will consist of an outer ring of 12, an inner ring of 6 and one electron at the centre. 

Following the above jprocedure we can construct a table of Mendeljev and compare it with the results of the 'exact' 
numerical simulations.EJ This is done by finding the distribution of the electrons when they are arranged in what we 
may consider to be the simplest way, i.e. when the number of rings is a minimum. The number of electrons in the 
outer ring Ni will then be determined by the equation 

N-N 1 =f(N 1 ). (2.9) 

The value of Ni as obtained from this equation is not an integer and consequently we have to take the integral part 
of this value. To obtain 7V2, the number of electrons in the second ring, we solve 

N-N 1 -N 2 = f(N 2 ). (2.10) 

We continue in this way until there remain less than 6 electrons. This procedure results into Table || This table of 
Mendeljev for the Thomson model is compared with the 'exact' table of Ref. [|. The configurations for which the 
Thomson model gives the wrong results are printed in italics. Notice that this model is capable to predict most of 
the configurations correctly. For systems with many electrons (N > 35) the Thomson model starts to fail. However 
even in this case the number of rings is still predicted correctly until N < 50. The reason for this difference is that 
in the Thomson model the configurations are found by using a stability argument, while in Ref. [l] they were found 
by a minimalization of the energy using Monte Carlo simulations. Here we will not compare the energy as obtained 
from the Thomson model with the result given in Ref. |. We found that the energy in the Thomson model for TV > 6 
deviates strongly from the 'exact' results which is due to the assumption that all inner electrons are placed at the 
centre. 



III. EXTENSIONS OF THE THOMSON MODEL 
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A. Effect of the confinement potential 



It turned out that for parabolic confinement and Coulomb repulsion the maximum number of particles on the 
inner ring is 5. Now we want to investigate how the confinement potential influences the possible configurations. For 
simplicity, we take n' = 1 for the inter-particle repulsion, i.e. Coulomb repulsion. 



From Eq. (2.7) we obtain the minimal value for p which stabilizes a ring of N out electrons 



P> Y (L - h) - —S N + T (2 - „)tf + - (3.1) 

and this inequality must hold for every k. Notice that only for n = 2, i.e. parabolic confinement, this stability 
condition does not depend on the angular velocity £1. For a linear confinement potential (n = 1) each configuration 
will become unstable for sufficient large f2, while for n > 2 every configuration can be stabilized if the system rotates 
fast enough. To get an idea of the influence of the confinement potential on the configurations we limit ourselves to 
the case SI = 0. Following the same procedure as used in previous section we find the table of Mendeljev (Table ||). 
The results presented in Table |l| are a good first approximation to the 'exact' results for the considered iV-values. 
Figure [| shows the maximum number of electrons on the inner, second, sixth ring as function of n. Notice that 
for a linear confinement potential the inner ring can support a maximum of 3 electrons which is much smaller than 
in the case of a parabolic confinement. For confinement potentials with n > 2 it's just the opposite, more electrons 
can be fitted on the rings without destabilizing them. 

In the limit of n — > oo the confinement potential becomes a hard wall potential, i.e. V(r) — (r < A), oo (r > A). 
We found numerically that the maximum number of particles on each of the rings keeps increasing with n. This 
clearly signals the breakdown of the Thomson model in this limit because from 'exact' numerical simulations (see 
Ref. |l|) with a hard wall circular potential it is known that, with increasing number of electrons, inner rings are 
formed. We can understand this as follows. In the Thomson model the groundstate configurations are found using a 
stability argument, so only inner rings are formed to stabilize the outer electrons. For a hard wall potential, there is 
no confining force for r < A, thus no inner electrons are needed. All the electrons form one ring at r — A, but this is 
not the configuration with the minimum value of the energy. „ 

For completeness, we mention that recently Farias and PeetersQ studied a system of 4 electrons in a Coulomb type 
of confinement and found different configurations depending on the strength of the Coulomb confinement potential. 
This together with the present results indicate that the number of electrons on each ring are not universal but can 
be strongly influenced by the type and strength of the confinement potential. 



B. Effect of the inter-particle interaction potential 

Now we investigate the effect of the functional form of the interaction potential between the electrons on the 
groundstate configurations. For simplicity, we consider here a parabolic confinement potential (n = 2) 



First we can extend the inter-particle potential to a logarithmic interaction. To obtain Eqs. (2.4) and (2.7) we 
used 1 [r n as the inter-particle potential. These equations are obtained from the law of Newton which contains the 
force. Because the derivative of In r equals 1 /r these equations are therefore also valid for a logarithmic inter-particle 

and 1(0) = 1. 



In Table III, we give the table of Mendeljev as obtained from the generalized Thomson model using a logarithmic 



and a 1/r interaction, and compare the results with the one for a Coulomb interaction (Table III). For the 1/r 
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interaction we find that with increasing number of electrons more electrons have to be placed at the centre than for 
the case of Coulomb interaction in order to stabilize the outer ring. Therefore the resemblance with the 'exact' table 
of Mendeljev will not be as good as for the case of Coulomb interaction. Indeed, the effect of the assumption that all 
inner electrons are placed at the centre will be larger. For the logarithmic interaction the opposite is found. Fewer 
electrons have to be placed at the centre and in this case the Thomson model should work very well. Figure [| shows 
the maximum number of electrons on the inner, second, sixth ring as function of n' . Notice that the maximum 
occupation number of the different rings decreases with increasing n' . 

For n! — » oo we found numerically Ni = 4 and for all the outer rings N — 6. Also in this limit we do not expect that 
the Thomson model gives the correct results, because the inter-particle interaction becomes extremely short range, 
i.e. they are delta function like. Classically the electrons can sit very close to each other and thus will occupy the 
r ~ region where the confinement potential is zero. As a consequence, such a system is similar to the case of an 
infinite 2D Wigner lattice where the particles form a hexagonal -lattice, i.e. each electron has 6 neighbours, and no 
ring structure is expected. Similar results were recently obtainedu for a screened Coulomb interaction e~ Kr /r using a 
molecular dynamics simulation approach. With increased screening, i.e. larger k, structural transitions were found in 
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which the ring configurations changed abruptly and in the limit of very large k a 2D Wigner type lattice was formed 
in the centre of the classical 2D atom. 



IV. EIGENFREQUENCIES INDEPENDENT OF THE NUMBER OF ELECTRONS 

Figure |l| suggests that in the case of parabolic confinement and Coulomb repulsion there are three eigenfrequencies 
independent of the number of electrons (see also Ref. |||). The existence and value of these eigenfrequencies can 
be obtained analytically. 1) For any axial symmetric system the system as a whole can rotate, which leads to an 

eigenfrequency u> — 0. 2) The Hamilton equation of motion yields v X i = —nx; L [x 2 + y 2 ) n ^ 2 1 +n' Y^j-^ ( x i~ x j) l T % +2 ■ 

Now consider the centre of mass R = ^ i ft which satisfies the differential equation 

d 2 R x _ . _ sr^ / 2 , 2N™/ 2 - 1 



dt 2 4^ 



and of course the same for R y . Notice that Eq. (4.1) is independent of n'. Thus only for a parabolic confinement 
potential the above equation reduces to d 2 R x /dt 2 — —2R X and a twofold degenerate vibration of the centre of mass 
with eigenfrequency w = v2 is obtained. This frequency is independent of the number of electrons and independent 
of the inter-particle potential which is a consequence of the generalized Kohn theorem.EI 3) For the mean square radius 
R 2 =12i{x 2 + yi) we find 

^ = -(n + n') £ (x 2 + y 2 )" /2 + n'H + 2T, (4.2) 

i 

with T = ^\ (xf + yf\ the total kinetic energy. For parabolic confinement, i.e. n = 2, there is a breathing mode with 
frequency u = \J1 + n 1 which is independent of the number of electrons. The existence of the breathing mode does 
not depend on the functional form of the inter-particle potential, but its value does. 



V. CONCLUSIONS AND SUMMARY 



The Thomson model, a model for classical 2D atoms, was investigated and the groundstate configurations were 
obtained. The results are valid for repelling particles with arbitrary mass and charge. First, we deduced within 
the Thomson model a table of Mendeljev for parabolic confinement aad Coulomb repulsion between the electrons 
and compared it with the results from 'exact' Monte Carlo simulations.!!! The Thomson model correctly predicts the 
configurations for N < 36. This table was constructed using stability arguments: the eigenfrequencies of the 2D 
atom are calculated and determined when one of them becomes imaginary in order to find the maximum number of 
electrons on the rings. 

Knowing that the Thomson model is a rather good model to predict the configurations of classical 2D atoms we 
investigated the effect of a r"-confinement potential. Only for parabolic confinement the configurations do not depend 
on the angular velocity for rotation of the total system. For a linear confinement potential the configuration becomes 
unstable if the system rotates too fast, while for n > 2 every configuration may be stable if the system rotates fast 
enough. We also found that the maximum number of allowed electrons on each ring is an increasing function of n, 
i.e. it increases with the steepness of the confinement potential. 

We extended the Thomson model to a 1 /r n -type and to logarithmic interaction between the electrons and found the 
groundstate configurations. Here we found that with increasing n', i.e. for shorter range of inter-particle interaction, 
fewer electrons can be put on a ring. 

We showed that there are three eigenfrequencies which are independent of the number of electrons in the case of 
parabolic confinement. The zero frequency mode which is a consequence of the axial symmetry of the system. We 
found that only for a parabolic confinement potential the frequencies of the vibration of the centre of mass and the 
breathing mode are independent of the number of electrons. The value of the latter does depend on the functional 
form of the inter-particle potential. 
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FIG. 1. The eigenfrequencies squared for N — 2,3, 4, 5, 6 for Q — 0. 
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FIG. 2. The number of electrons f(N out ) needed in the centre of the system in order to stabilize a ring of N out electrons. 
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FIG. 3. The maximum number of electrons on the inner, second, sixth ring as function of the power of the confinement 
potential (n) for the case of Coulomb repulsion between the electrons. 
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FIG. 4. The maximum number of electrons on the inner, second, sixth ring as function of the power of the inter-particle 
interaction (n') for the case of a parabolic confinement potential, n = corresponds to a logarithmic interaction. 



TABLE I. Table of Mendeljev for classical 2D atoms. The results as obtained from the Thomson model are compared with 
the 'exact' results from the Monte Carlo simulations of Ref. ll 



N 


Thomson 


Monte Carlo 


N 


Thomson 


Monte Carlo 


1 


1 


1 


26 


3,9,14 


3,9,14 


2 


2 


2 


27 


4,9,14 


4,9,14 


3 


3 


3 


28 


4,10,14 


4,10,14 


4 


4 


4 


29 


5,10,14 


5,10,14 


5 


5 


5 


30 


5,10,15 


5,10,15 


6 


1,5 


1,5 


31 


5,11,15 


5,11,15 


7 


1.6 


1,6 


32 


1,5,11,15 


1,5,11,15 


8 


1,7 


1,7 


33 


1,6,11,15 


1,6,11,15 


9 


2,7 


2,7 


34 


1,6,12,15 


1,6,12,15 


10 


2,8 


2,8 


35 


1,6,12,16 


1,6,12,16 


11 


2,9 


3,8 


36 


1,7,12,16 


1,6,12,17 


12 


3.9 


3,9 


37 


2,7,12,16 


1,7,12,17 


13 


4,9 


4,9 


38 


2,7,13,16 


1,7,13,17 


14 


4,10 


4,10 


39 


2,8,13,16 


2,7,13,17 


15 


5,10 


5,10 


40 


2,8,13,17 


2,8,13,17 


16 


5,11 


1,5,10 


41 


2,9,13,17 


2,8,14,17 


17 


1,5,11 


1,6,10 


42 


3,9,13,17 


3,8,14,17 


18 


1,6,11 


1,6,11 


43 


3,9,14,17 


3,9,14,17 


19 


1,6,12 


1,6,12 


44 


4,9,14,17 


3,9,14,18 


20 


1,7,12 


1,7,12 


45 


4,9,14,18 


3,9,15,18 


21 


2,7,12 


1,7,13 


46 


4,10,14,18 


4,9,15,18 


22 


2,7,13 


2,8,12 


47 


5,10,14,18 


4,10,15,18 


23 


2,8,13 


2,8,13 


48 


5,10,15,18 


4,10,15,19 


24 


2,9,13 


3,8,13 


49 


5,11,15,18 


4,10,15,20 


25 


3,9,13 


3,9,13 


50 


1,5,11,15,18 


4,10,16,20 
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TABLE II. Table of Mendeljev for the Thomson model for various confinement potentials r™. 



N 


n=l 


n=2 


n=3 


n=10 


N 


n=l 


n=2 


n=3 


n=10 


3 


3 


3 


3 


3 


14 


1,4,9 


4,10 


4,10 


14 


4 


1,3 


4 


4 


4 


15 


1,4,10 


5,10 


4,11 


15 


5 


1,4 


5 


5 


5 


16 


1,5,10 


5,11 


5,11 


16 


6 


1,5 


1,5 


6 


6 


17 


1,5,11 


1,5,11 


5,12 


17 


7 


1,6 


1,6 


7 


7 


18 


1,6,11 


1,6,11 


6,12 


1,17 


Q 


2,b 


1,7 


1,7 


Q 

o 


1 CI 

ty 


o C 11 

z,o,±± 


1 c n 
1,0, Iz 


7 10 
1 ,lZ 


2,17 


9 


2,7 


2,7 


1,8 


9 


20 


2,7,11 


1,7,12 


7,13 


3,17 


10 


2,8 


2,8 


1,9 


10 


21 


2,7,12 


2,7,12 


1,7,13 


4,17 


11 


3,8 


2,9 


2,9 


11 


22 


2,8,12 


2,7,13 


1,8,13 


4,18 


12 


3,9 


3,9 


2,10 


12 


23 


3,8,12 


2,8,13 


1,9,13 


5,18 


13 


1,3,9 


4,9 


3,10 


13 


24 


3,8,13 


2,9,13 


1,9,14 


6,18 



TABLE III. The groundstate configurations for logarithmic, Coulomb and 1/r 2 interaction between the particles, as obtained 
within the Thomson model, as a function of the total number of particles N. 



N 


— In r 


1/r 


1/r 2 


N 


— lnr 


1/r 


1/r 2 


2 


2 


2 


2 


27 


2,9,16 


4,9,14 


5,10,12 


3 


3 


3 


3 


28 


2,9,17 


4,10,14 


5,10,13 


4 


4 


4 


4 


29 


2,10,17 


5,10,14 


1,5,10,13 


5 


5 


5 


5 


30 


2,10,18 


5,10,15 


1,6,10,13 


6 


6 


1,5 


1,5 


31 


3,10,18 


5,11,15 


1,6,11,13 


7 


1,6 


1,6 


1,6 


32 


3,11,18 


1,5,11,15 


1,7,11,13 


8 


1,7 


1,7 


1,7 


33 


3,11,19 


1,6,11,15 


2,7,11,13 


9 


1,8 


2,7 


2,7 


34 


4,11,19 


1,6,12,15 


2,8,11,13 


10 


2,8 


2,8 


2,8 


35 


4,12,19 


1,6,12,16 


2,8,11,14 


11 


2,9 


2,9 


3,8 


36 


4,12,20 


1,7,12,16 


3,8,11,14 


12 


2,10 


3,9 


3,9 


37 


5,12,20 


2,7,12,16 


3,8,12,14 


13 


3,10 


4,9 


4,9 


38 


5,13,20 


2,7,13,16 


3,9,12,14 


14 


3,11 


4,10 


5,9 


39 


5,13,21 


2,8,13,16 


4,9,12,14 


15 


4,11 


5,10 


5,10 


40 


6,13,21 


2,8,13,17 


5,9,12,14 


16 


4,12 


5,11 


1,5,10 


41 


6,14,21 


2,9,13,17 


5,10,12,14 


17 


5,12 


1,5,11 


1,6,10 


42 


6,14,22 


3,9,13,17 


5,10,13,14 


18 


5,13 


1,6,11 


1,6,11 


43 


1,6,14,22 


3,9,14,17 


1,5,10,13,14 


19 


6,13 


1,6,12 


1,7,11 


44 


1,7,14,22 


4,9,14,17 


1,5,10,13,15 


20 


6,14 


1,7,12 


2,7,11 


45 


1,7,15,22 


4,9,14,18 


1,6,10,13,15 


21 


1,6,14 


2,7,12 


2,8,11 


46 


1,7,15,23 


4,10,14,18 


1,6,11,13,15 


22 


1,7,14 


2,7,13 


3,8,11 


47 


1,8,15,23 


5,10,14,18 


1,7,11,13,15 


23 


1,7,15 


2,8,13 


3,8,12 


48 


1,8,16,23 


5,10,15,18 


2,7,11,13,15 


24 


1,8,15 


2,9,13 


3,9,12 


49 


1,8,16,24 


5,11,15,18 


2,8,11,13,15 


25 


1,8,16 


3,9,13 


4,9,12 


50 


2,8,16,24 


1,5,11,15,18 


2,8,11,14,15 


26 


2,8,16 


3,9,14 


5,9,12 


51 


2,9,16,24 


1,6,11,15,18 


3,8,11,14,15 
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